Bogoliubov transformation for distinguishable 

particles 



M. Holzmann and F. Laloe 
Laboratoire de Physique de FENS, 
LKB et LPS, 24 rue Lhomond, F- 75005 Paris, France 

February 1, 2008 
Abstract 

The Bogoliubov transformation is generally derived in the context 
of identical bosons with the use of "second quantized" a and cd opera- 
tors (or, equivalently, in field theory). Here, we show that the transfor- 
mation, together with its characteristic energy spectrum, can also be 
derived within the Hilbert space of distinguishable particles, obeying 
Boltzmann statistics; in this derivation, ordinary dyadic operators play 
the role usually played by the a and cd operators; therefore, breaking 
the symmetry of particle conservation is not necessary. 

The Bogoliubov transformation Q 0] is an essential tool in the theory of 
Bose-Einstein condensation of identical bosons^]. It modifies the quadratic 
energy spectrum of free particles into a quasiparticle spectrum which in- 
cludes a linear variation for small momenta; this feature is generally asso- 
ciated with the existence of phonons and, since it introduces a non zero 
minimum value for the ratio between the energy and momentum, it allows 
a natural derivation of the notion of critical velocity (a maximum veloc- 
ity for the system to remain superfluid). Usually, the mathematics of the 
Bogoliubov transformation is performed within the formalism of creation 
and annihilation operators (often called "second quantization" for historical 
reasons); assuming that the system is entirely made of identical particles, 
one then uses a formalism which automatically ensures a full symmetriza- 
tion of the state vector. Nevertheless, the notion of phonons has a much 
broader scope in physics than just identical quantum particles; it is even 
often discussed in the context of classical systems, solids or even fluids. One 
can therefore wonder whether it is possible to re-derive the Bogoliubov spec- 
trum in a context where the particles are considered as distinguishable and 
where, as a consequence, the effect of exchange operators remains completely 

1 Historically, the first introduction of the mathematical transformation seems to be the 
work of Holstein and Primakoff in 1940 S , in the context of magnetic systems. 



explicit. The purpose of the present article is to show how this is indeed pos- 
sible. Another motivation for a mathematical derivation of the Bogoliubov 
transformation within the Hilbert space of distinguishable particles arises 
for the use of Ursell operators in statistical mechanics [||, a formalism in 
which the symmetrization of the states is not introduced implicitly from the 
beginning of the calculations, but explicitly and at a later stage with the 
help of exchange cycles. 

1 Hamiltonian 

The hamiltonian of the problem is: 

where N is the number of particles, m their mass, Pj the momentum of par- 
ticle numbered i, and V(i,j) the interaction energy of particles numbered i 
and j. The simplest assumption is to take the matrix elements of this inter- 
action potential as constant, provided they satisfy momentum conservation 
(otherwise they of course vanish) : 

< i : k , j : k | V(i,j) \ i : k + q , j : k' - q >= g (2) 

where g is the coupling constant, inversely proportional to the volume of the 
systemF). The hamiltonian can then be written: 

^ = E^T + |E E |<:k, j:k'>«:k + q> j:k'-q| (3) 
i=1 #3 k,k',q 

In this expression the interaction term contains, first, the forward scat- 
tering terms q = which can be written: 

f EE i i : k y<i : k i 1 j ' : ^><i '■ k 'i 

i#i k,k' m 

=!EEi* :k ><* :k i^ :k '><J :k 'i-fiv 

iJ k,k' 

or simply: 

fiV(JV-l) (5) 

2 In the thermodynamic limit, g itself tends to zero, but products such as gN keep a 
finite value. 
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As for the q 7^ terms, they never contain 4, or even 3, vanishing momenta; 
the terms containing 2 vanishing momenta can include them, either in the 
same side of the operator, or in opposite sides: 



Yl [l i : ' 3 ■ >< * : k . 3 ■ -k I +h.c. 

I [1 * : ' j : k >< * : k > j : 1 +h - c - 



(6) 



where /i.e. is for Hermitian conjugate. It is convenient to express the second 
line0 of this expression as a function of simpler (diagonal) operators by re- 
writing it in the form: 



9 12Y1 [l * : k ' J -Oxi-.k, j : 



+ W a 



(7) 



where W as . is the (antisymmetrical) interaction operator: 

W as . : 



gYsYl i 1 - P ™chXi,j)) I i ■ k , j : >< i : k , j : | (8) 

i^j k 



Here, Pexch.(hj) 1S the exchange operator of particles i and j (the condition 
k 7^ in the summation can be released since the corresponding term van- 
ishes); another equivalent expression of W as _ can be obtained by applying 
the exchange operator to the rightQ 



W as . = -9 Y Yl I 1 : ' j : k >< * : ' j : k 

i+j k 



exch. [}: J , 



(9) 



(in passing, we see that W as . is Hermitian). Now, the first operator in (^) 
can be simplified into: 



'EE [1- 



k >< i : k I (g> I j : >< j : 



i,j k^o 



(10) 



(the constraint % 7^ j can be released since every term % = j in the summation 
is zero, due to the orthogonality of single particle states), or again: 



gN Y n k = gN N e 

k^O 



(11) 



For this operator in the second line, interchanging the dummy indices i and j is 
equivalent to an hermitian conjugate operation; for this term, we can therefore ignore the 
h.c. and just replace g/2 by g . 

4 In (Q), (||) or (^), interchanging the dummy indices i and j is equivalent to interchang- 
ing the states labelled by and k. 
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where is the population of the excited state labelled by momentum k: 

n k = ^ | i : k >< i : k | (12) 

i 

Nq the population of the ground state (for the ground state, we use a capital 
letter to emphasize that it has an extensive population, but iVo is in fact the 
same operator as no): 

N = I * : >< * : I ( 13 ) 

i 

N e as the operator associated with the total number of excited particles: 

N e = J2^ = N -No (14) 
Finally, the interaction terms with 2 vanishing momenta can be written: 



^^[|t:0,j:0><j:k,i : -k|+/i.c.] +gN N e + W as . (15) 



Last, the only interaction term that we have not yet included corresponds 
to the interaction between particles in excited states; we call this operator 

Vee- 



Vee = ^12Yl E [M : k,j : k' >< i : k + q,j : k'-q| +/i.c. 



ift k,kV0 q^O,-k,+k' 



(16) 



but, in what follows, we will merely neglect its effect; this is because the 
system will be supposed to be at sufficiently low temperature and to be 
sufficiently dilute so that most of the particles remain in the ground state; 
interaction effects proportional to the square of the excited state populations 
are then negligible. 

To summarize, we have obtained the following expression for the Hamil- 
tonian: 



H =^(e fe +^o)n k + |AT(AT-l) 

k^O 

+| EE [i 1 ■ °> j : >< i ■ k >j ■ - k i +/i - c -i (17) 

ijij k^O 
+ V ee + W as . 

where the free particle energy is defined, as usual, by: 

« = -ST (18) 
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2 Another expression of the hamiltonian 



We introduce in this section a new hamiltonian H which in a second step, 
we will identify with H term by term. 

2.1 Introducing new creation and annihilation operators 

We now define the operator by: 

Ak = ^^{«k I * : >< i : k | +/3 k | i : -k >< i : |} ( 19 ) 



N 

i 

where we assume that a k and /3k are real numbers, which are for the moment 
not fixed, but which are supposed to be even functions of the vector k: 

«k = a-k ; /3 k = P-k (20) 

j4 k is defined as an operator which removes a momentum hk from the system 
by, either transferring one particle from momentum k to zero, or from zero 
to — k; it is a single particle operator which, if restricted within the totally 
symmetric part Ss of the Hilbert space, could be expressed (through the 
well known expression of single particle operators) as: 

^k = {"kao a k + /3kaoal k | (21) 



N 

with the usual notation for the creation and annihilation operators a k and 

4- 

The Hermitian conjugate of is equal to: 

A l = TTfi Yl W I j : k >< j : I +/3k I j : >< j : " k 1^ (22) 

j 

We can now calculate the commutator of A^ and A^ k term by term. For 
instance, the commutator of | i : >< i : k | with | j : k >< j : | is zero 
except if i = j; assuming that this is the case, the commutator is given by: 

| i : >< i : k | i : k >< i : | - | i : k >< i : | i : >< i : k | (23) 

which is merely the difference! i : >< i : | — | i : k >< i : k |. Taking 
the sum of all four terms in a similar way gives (the crossed term in x /?k 
vanish) : 



= N' 1 ^ {«k [l i : >< i : I - I i : k ><« : k 

i 

+ Pi [j i : -k >< i : -k | - \i : >< i : 



(24) 
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or: 



or again: 



^k , K 



N ° k iV 



(«k - AJ + A< — ^ a k 



(25) 



(26) 



We now calculate the product 4^ ^ or an y k ^ ^> ^ ne terms i = j 
provide: 



iV" 1 ^2 \al | i : k >< % : k | | z : >< i : 



while the terms i ^ j contain a summation of the expression: 
iV _1 |«k I * : 0) 3 ■ k >< i : k, j : | | « : -k, j : >< i : 0, j : -k 



(27) 



+"k x /? k 



j : 0, j : >< i : k, j : — zk | +/i.c. 



} 



(28) 

where, as above, h.c. is for Hermitian conjug atefl We can now use the same 
method as in § to distinguish, within the first line of (|28|), a "diagonal" 
part and an antisymmetrical part B as _(k): 



— [a£ N n k + fa N n_ k ] + B as . (k) 



(29) 



with the definition: 
1 



B as .(k) = -- J2 [1 - PexchXhO)} { [«k I * = k, j : X i : k, j : 

+ ^|i:0,i:-k ><i:0j':-k |]} 



or, equivalently, by applying the exchange operator on the other side: 
# as .(k) = -— { [ a l I * : °. 3 ■ k >< i : 0, j : k | 



+ Pi | » : -k, j : >< i : -k, j : 



| [1 - Peach. (*>.?')] 



(30) 



(31) 



Finally, we obtain: 

4^ = 4 nk(^)+/«^(l + n_ k ) 



+N 1 a k /? k < 



^ | j : 0, j : >< i : k, j : -k | +h.c. \ + P as .(k) 



(32) 



3 For this hermitian conjugate term, we have interchanged the dummy indices i and j. 
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2.2 Identification of two expressions 

Les us introduce a new hamiltonian H' by: 

H' = J2^k + XnN + A e iV e (33) 

where a^, /3k, e/T, and A e are for the moment free parameters - in a 
second step they will be chosen in order to make the new hamiltonian equal 
to the initial Hamiltonian: 



H =^(e fe + 5^o)n k + |iV(iV-l) 



o^v- ( 34 ) 
+ 2 ^ ^ [l * : °' j 1 ° * : k ' j : " k 1 

In addition, and will have commutation relation that are similar to 
those of ordinary creation and annihilation operators, 
(i) first condition (commutation relation); if: 

.2 p2 



ai - ft = 1 (35) 



(for any value of k), relation (|2q ) shows that the commutator of and 
A^ is equal to one in the limit of very low temperatures and very dilute 
systems (when almost all the particles are in the ground state, N e -C N). 
Relation ( |35|) is automatically fulfilled with the following choice of the two 
parameters ctk and /3k as a function of a single parameter 

a k = cosh £ k / 3g N 
/3 k = sinh £ k 

(ii) second condition (identification of the main interaction terms) ; from 
(0), © and (jH) we get that the terms in | i : 0, j : >< i : k, j : -k | 
can be made identical if we set: 

e^OkAc = ^— (37) 
which, through ([?(]), is equivalent to: 

^sinh2^ k = ff iV (38) 

(iii) third condition (kinetic terms in n k ); the identification of the terms 
which, in (fj^), are linear in the excited population operators n k (or n_ k ) 
provides the condition: 

r 2 i + /v . , 2 iv \ 

e fc i "k jy + (A-k) f = e fc + 5^0 (39) 
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or, through the relation N = Nq + N e : 



2 1 + N-N e , 2 N-N e 
efc i « k r T + (A-k) 



A r 



N 



e k + gN- gN e 



(40) 



In this equation, N is a number while N e is an operator; term by term 
identification then provides the two conditions: 



Gk 



or 



;i + iv- 1 ) + (/3_ k r 



e k + gN 



Assuming that AT 3> 1, and taking into account the parity relation 



well as definition (36), we can write the former in the form : 



efc («k + £k) = e k cosh 2£ k = e k + gN 



(41) 



MS 



(42) 



an equality which, together with (|38|) , will provide the Bogoliubov quasipar- 
ticle spectrum. In the second line of (]4l|), we have intentionally left in both 
sides the operator N e ; in this way we emphasize that, in the expression of 
H , this term appears as a product of N e by the population operator n k of 
the excited state k, in other words as a second order correction in N e /N 
which can be neglected in the limit of low temperatures and very dilute 
systems. Therefore, the major constraint of the identification is contained 
in ( f42| ) and, from now on, we will leave aside the second condition of (fill). 

(iv) terms in A" and N e ; in d32"D, we have not yet included the effect of 
the term (/3 k ) 2 Nq/N which, when Nq is replaced by A^ — N e , provides: 



N(N - 1) = J2 % (A<) 2 (i " ^) + X nN + 



X,N P . 



(43) 



Term by term identification then provides^ 



(44) 



and: 



(45) 



We will neglect A e in what follows. 

6 Since, according to (|47|), /3k will be proportional to l/e^ when k is large, the 



equation ( |43[ ) and following are linearly divergent at infinity. We discuss in § 3.1.2 
this divergence can be eliminated 



sums m 
how 
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(v) antisymmetric terms B as . and W as ; the initial hamiltonian contains 
the operator W as . while H contains the operator: 

Was. = ^2^B asik) (46) 

k 

where B as ^ is defined in (30) or (|3l|) . We note that both these terms con- 
tain two particle antisymmetrizers [1 — P eX ch.hj)] an d will therefore always 
vanish when multiplied (on any side) by the N particle symmetrization op- 
erator Sn; their contribution is therefore exactly zero if the particles in 
the system are identical bosons. By the same token, the same remains true 
if the particles are distinguishable but in their ground state, which is also 
completely symmetrical (it is actually exactly the same as for bosons). The 
property obviously extends to any excited state having the same permuta- 
tion symmetry, but not necessarily for energy states which correspond to 
other representations of the permutation group; in general, there is no rea- 
son why the difference W as . — W as should play no role for distinguishable 
particles^. 



2.3 Bogoliubov spectrum 

From (^8|) and ( |42|) we get by taking the ratio: 



tanh2g k = y 47 

e k +gN 

which fixes the value of to determine we express cosh 2£ k as a function 
of this result: 

= 1 - tanh 2 2^ k = 6k {6k + (48) 
cosh 2 2& (e k +gN) 2 

which, combined with ((42|), gives: 

^ = Vefc (ejfc + 2gN) (49) 

This is the well known Bogoliubov result for the energy of the quasi-particles; 
actually relations (^7|) and (49) are exactly the basic relations obtained in 
the usual calculation in terms of annihilation and creation operators. 



7 Using the parity of /3k and interchanging the dummy indices i and j in (]3CJ) allows to 
show that W as . is equal to the sum over k of the product e~k (qJ + /?£) by the same sum 
over i ^ j which appears in (^); the strict equality of W as . and Was. would require the 
condition JZk^ ( a k + fii) = gN, which is reminiscent of the second condition (fid]), but 
introduces convergence problems. Requesting a strict equality of the operators H and H 
in the whole Hilbert state of Boltzmann particles would then lead to difficulties. 
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With the above relations, the two Hamiltonians H and H can be iden- 
tified, with three approximations: 

* (Al) we assume that the difference W as , — W as can be ignored, which 
is exact for totally symmetric states. 

* (A2) we ignore the effect of V ee , corresponding to the interactions 
between excited particles. 

* (A3) we ignore the condition expressed in the second line of (|4l|). 
The third approximation is consistent with the second and remains valid 

at low density and temperature, exactly as in the usual derivation for bosons. 
The first is of different nature since it is necessary only for Boltzmann parti- 
cles; we have already mentioned that the operator W as . — W as plays no role 
for the ground state of the system, as well as for all states which are fully 
symmetric - for instance for all states that are obtained from the ground 
state by action of any product of creation operators A^s, since these oper- 
ators do not change the permutation symmetry of the states. 



3 Discussion 

We now discuss more precisely how the operators and can be used 
to construct eigenstates of the hamiltonian. 



3.1 Commutation relations 



One of the commutation relations between the new creation and annihilation 
operators has already been studied above, and led us to an approximation 
which is consistent with approximation (A2): 

* (A4) we ignore, in (26), the terms in N e /N (as well as those in n^/N, 
but the latter tend to zero in the thermodynamic limit and therefore raise 
no question). 

But we also have to study the situation for different values of k as well 
as other commutations relations. We see in the definition ( |l9|) of that 
its commutator with A / will, first, contain the commutator: 



i : >< i : k I, I i : >< i : k 



(50) 



which vanishes since both products of the operators do so (orthogonality of 
single particle states if k and k are both different from zero). Similarly: 



i : -k >< i : |, | i : -k >< i : 
Two other commutators do not vanish, for instance: 
: -k >< i : |, | i : -k' >< i : | 







i : 



-k >< i : k 



>< i : 



(51) 



(52) 
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so that we obtain the result: 



[Ak, A k '] = <5 k _ k / (a k /3_ k - a-k /3k) ^ 

-^VlK-k'xKkl+^Vl.-kXrk'l (53) 

i i 

The parity relation ( p0| ) ensures that the main term in the right hand side, 
proportional to Nq/N, vanishes exactly. The remaining terms, on the second 
line, are as V ee "excited-excited terms" which act only on excited particles, 
and leave them excited; it is therefore natural to introduce one more as- 
sumption: 

* (A5) the operators in the second line of d53| ) can be neglected in our 
calculation, so that the A^s get the usual commutation relations of annihi- 
lation operators for orthogonal states. 

We do not have to study the commutation relations of A^ and Aj, , since 
they can be obtained from (|53|) by Hermitian conjugation, but we have to 
study the commutator of A^ and A^, . The calculation from (19) and ( p2[ ) is 



actually very similar to that which leads to ( |53| ) and provides: 

-^y|,:k'>< i: k| + Mr| j: _ k><l: _ k '| ^ 

i i 

When k = k , this relation has already been studied; when k ^ k', we 
introduce the additional assumption, similar to (A5): 

* (A6) the operators in the second line of fl54| ) can be neglected in our 
calculation since they also correspond to "excited-excited terms" . 

3.1.1 effect of the operators A^ and A^ 

For completeness, we recall here a last approximation, which has already 
been discussed above and which is consistent with all preceding approxima- 
tions (valid if the gas is very dilute and at a very low temperature): 



* (A7) we ignore, in the expression (33), the effect of the term in N e 
(the term in N creates no problem since this number is unchanged under 
the action of all operators introduced; this term is studied below). 

Assume then that | &e > is an eigenstate of the Hamiltonian H, with 
eigenvalue E; we introduce the new ket | <&^(k) > as: 

| $^(k) >= A k | <S> E > (55) 

Now, since all usual commutation relations are satisfied by the opera- 
tors j4k, we can easily show that | ^^(k) > is another eigenstate of the 
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Hamiltonian H, with eigenvalue E — A^ then plays the role of a "lad- 
der" operator which changes the energy step by step, by decreasing values; 
similarly, will increase the energy eigenvalues. 

Let | $o > be ground state of the hamiltonian H. By action of A^ onto 
this ket we obtain, either a ket which is zero, or a ket in which the average of 
H has decreased by e*T , plus possibly some corrections related to the terms 
in the commutators that have been neglected. As noted above, | $o > is a 
common ground state to bosons and Boltzmann particles, since it is fully 
symmetric; so is A^ | $q >, which shows that the effect of W as . — W as 
on this ket is strictly zero. We then just have to deal with approximations 
(A2) to (A7), which amount to assuming that the ratio N e /N has a small 
average value in the ground state. If the gas is sufficiently dilute, all the 
corrections in N e /N will not be able to make up for the decrease in energy 
e/T and change it into an increase of energy. But no state has an average 
energy below that of the ground state, so that one necessarily has: 



This equation is valid for any value of k, except very small values for which 
the density corrections in N e /N may be comparable to the decrease in energy 
e^T, so that the cancellation of the ket is no longer a necessity. We have 
obtained in this way a set of equations that the ground state of the system 
has to obey. 

Now, applying any power of the operator A^ to this ground state will not 
change the permutation symmetry of the state and therefore allow us to still 
ignore the effect of W as . — W as _ ; we will therefore obtain good approximations 
to energy eigenstates, as long as not too many particles are excited in the 
operation. The Bogoliubov energies therefore play the role of quasiparticle 
energies, as in the usual derivation for bosons, which is not surprising since 
the completely symmetric subspace £s of the Hilbert space of Boltzmann 
particles remains invariant under all the operators considered. 

3.1.2 ground state energy 

Finally, it is interesting to come back to the term in \n that we have obtained 
in (H|). We have: 



A k | $ >= 



(56) 




(57) 



k 



k 



where, if we use (45), we can insert: 




(58) 
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which gives: 



X N = |(JV - 1) + + 2 9 N ) -e k -gN\ (59) 



When k tends to infinity, the expression under the sum is equivalent to: 

(gN? 



(60) 



A cancellation of first order terms occurs, but not of second order terms, so 
that the sum over all values of k remain (linearly) divergent. 

This divergence is well known and takes place in all derivations of the 
Bogoliubov transformation. A classical way to solve the problem || is to 
replace the coupling constant g (which is directly the matrix element of 
the potential) by its second order expansion as a function of the scattering 
length a: 

iirah 2 { Aitah 2 ^ 1 ) 

* = ^( 1 + — + (61) 

where V is the volume. Expressing the mean interaction term gN(N — l)/2 
as a function of a instead of g then introduced a counterterm in ( |59D which 
merely amounts to adding inside the sum (gN) 2 /ek, which pushes the can- 
cellation of orders up to second order, makes the sum convergent, and pro- 
portional to the integral: 



J x 2 dx | y / x 2 (x 2 + 2) - x 2 - 1 + } 



(62) 



The value of this integral turns out to be V128/15 and, finally, one gets for 
the following well-known result for the energy of the ground state: 



X N N = jrN(N -!)<! + — f-y — rr f (63) 




Another method to eliminate the divergence is to use a pseudopotential^ 
which directly contains the scattering length a in its matrix elements, so 



that any renormalization of the coupling constant such as (61) becomes 
uncecessary. In our calculations, we have replaced all matrix elements of 
the potential by the same constant g, but it would be possible to make a 



Not to be confused with an ordinary delta function potential; a real pseudopotential 
contains, in addition, a r derivation operatorQ^] and leads to a scattering lenght which is 
indeed proportional to its matrix elements, while it turns out the cross section associated 
with an ordinary delta function potential vanishes. 
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more careful calculation with a correct expression of the matrix elements 
of a pseudopotential; as shown by Castin Q, the method also allows one 
to recover (63). Beliaev [H] has studied systematically how the Bogoliubov 



energy spectrum with the scattering length a as an interaction parameter 
is recovered from a resummation of diagrams, as well as corrections to the 
Bogoliubov theory. 

CONCLUSION 

The mathematics of the Bogoliubov transformation can be performed 
within the space of states 8b of distinguishable particles, and leads exactly 
to the same formulas than in its fully symmetrical subspace £s; in other 
words, the proof can be generalized to include quantum states which are 
not necessarily completely symmetric by exchange and, as one would natu- 
rally expect, renders explicit additional conditions of validity. Even within 
the completely symmetric space 8s-, the proof is not only a mathemati- 
cal curiosity since it does necessitate a symmetry breaking of the number 
of particles; it therefore reveals more precisely what is behind the usual 
approximations that is made by replacing the operators ao and aj by c- 
numbers, a brutal approximation whose effects are not necessarily easy to 
control quantitatively. Here one gets a precise view of the exact list of opera- 
tors which have been neglected (approximations A2 to A7) - see for instance 



the terms in N e /N in (26) or those appearing in the second lines of ( J53| ) 
and ([34|), which arise precisely from the fact that ao and aj have not been 
treated as numbers; if one wished to push the approximation beyond the 
usual derivation of the Bogoliubov transformation, the consideration of the 
exact expression of these terms would be useful. Needless to say, during the 
derivation, we have also obtained trivial terms such as the interaction term 
between excited particles, V ee ,, which is ignored exactly in the same way as 
in the traditional derivation: it should be no surprise that, for bosons as 
well as for Boltzmann particles, the ratio N e /N should remain small for the 
Bogoliubov spectrum to be established. 

In a sense, the most interesting correction term is the operator W as , — 
W as , even if it has no effect on the whole class of states that are common to 
bosons and Boltzmann particles (those which can be obtained by the action 
of operators onto the ground state). Indeed, we know the low energy 
spectrum of a system of distinguishable particles is richer than that of a 
system of bosons, and here this fact is reflected in our calculation by the 
presence of these antisymmetric operators)^]. It is well known that the major 
difference between a Boltzmann system and a boson system is not to be 



9 In all calculations based on the Ursell operator formalism for bosons, where a com- 
plete symmetrization is applied in a second step, the two operators Was. and W as . give 
contributions which exactly vanish, so that the presence of these antisymmetric operators 
will create no problem. 
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found in the condensate, which is basically described by exactly the same 
many-body wave function in both statistics; it is in fact contained in the 
excitations, which are much more numerous for Boltzmann particles since 
the numbering of the particles which are excited becomes relevant. For non- 
interacting particles, this just associates much more entropy to excitation 
processes than for identical particles - this is actually the reason behind 
the instability of the condensate for Boltzmann particles at any non-zero 
temperature, as opposed to a Bose system where it remains stable until the 
Bose Einstein temperature is reached.. As soon as interactions are included, 
energetic effects also occur and, in our approach, they are reflected by the 
presence of the operators W as . and W as . 

LKB and LPS are Unites associees au CNRS et a PUniversite Pierre et 
Marie Curie, UMR 8552 and 8551. 

References 

[1] N.N. Bogoliubov, J. Phys. USSR 11, 23 (1947). 

[2] see for instance R.K. Pathria, "Statistical mechanics", Pergamon Press 
(1972), §11.2 or A.L. Fetter and J.D. Walecka, "Quantum theory of 
many-particle systems", Mc Graw Hill (1971), §35. 

[3] T. Holstein and H. Primakoff, Phys. Rev. 58, 1098 (1940). 

[4] P. Griiter and F. Laloe, J. Physique 5, 181 and 1255 (1995); M. Holz- 
mann, P. Griiter and F. Laloe, Eur. Phys. J. B 10, 739 (1999). 

[5] K. Huang and C.N. Yang, Phys. Rev. 105, 767 (1956) 

[6] K. Huang, "Statistical mechanics", Wiley (1963), chap. 13. 

[7] Y. Castin, seminar given in 1995 at ENS and unpublished notes. 

[8] S.T. Beliaev, ZH. Eksp. Teor. Fiz. 34, 417 and 433 (1958); JETP 7, 
289 and 299 (1958). 



15 



